TITLE OF THE INVENTION 
DIGITAL MODEM 

BACKGROUND OF THE INVENTION 
FIELD OF THE INVENTION 

The present invention concerns a digital modem for 
modulating/demodulating 1 /-1 binary signal taking as diffusion signal a signal 
of which group delay is proportional to the frequency 
BRIEF SUMMARY OF THE INVENTION 
[Problems to be Solved by the Invention] 

The spectrum diffusion method transfers a signal by performing the 
diffusion modulation using noise-like diffusion symbol and by diffusing the 
signal energy over a band width larger far than the information band width. 

Consequently, the carrier band is enlarged to reduce the energy density, 
allowing to reduce interference to the other radio waves. 

The reception side performs the diffusion demodulation using the same 
diffusion symbol as the transmission side, reducing radio interference based on 
the diffused signal, and can receive only the desired wave with a high SNR 
among many existing radio waves. 

In general, the spectrum diffusion method uses artificial noise series as 
diffusion symbol, but requires several conditions such as having a peak sharp in 
auto-correlation characteristics, cross-correlation characteristics evenly small, 
and a large number of available kinds of symbol, or others. 

Therefore, transmitters and receivers for generating diffusion symbols 
satisfying these conditions become complicated inconveniently. 
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Besides, to obtain correlation output by the receiver, it is necessary to 
synchronize correctly the timing of the reception signal and the reference signal 
generated at the receiver side. 

When the communication starts, as the timing is not known by the 
5 receiver side, the receiver side was required to perform considerably 
complicated synchronization control including operation to search for the 
connecting timing, operation for holding the obtained timing, or others. 

On the other hand, the cross-correlation between a signal h[i] of which 
group delay is proportional to the frequency and a signal h'p] which is hp] right 
1 0 and left inversed on the time axis presents the following characteristics. 

According to the definition thereof, the cross-correlation (in short, auto- 
correlation) of hp] and h[i] and h'p] and h'[i] is equal to a linear convolution of hp] 
and h'p] becomes an impulse. 

While, the cross-correlation of h[i] and h'p] is equal to the linear 
15 convolution of hp] and hp], or h'p] and h'p] and becomes substantially 0. 

Now, the reason thereof will be described as follows. 

hp] and h'p] is a kind of sweep signal as shown in Fig, 21 and Fig. 22. 

Ordinate axis of the graph represents the frequency (unit: Hz), while 
abscissa axis the time (unit: sample). 
20 Here, in order to simplify the analysis, the sweep signals hp], h'p] are 

approximated by a cosine signal of which frequency varies gradually as shown 
in Fig. 23 and Fig. 24. 

The sampling frequency is set to 1000 Hz, and the cross-correlation is 
determined for every 20 samples in correspondence to the frequency variation. 
25 In the following discussion, the signal amplitude scaling is not considered 

particularly, in order to simplify the description. 
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In the determination of cross-correlation of h[i] and h[0 , the value 
thereof becomes 40 when the time difference tau is 0 as shown below. 
Zh[i] h[i-faa]^=Zh[i]h[i] 

i i = 0 
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20 On the other hand, when the time difference tau is 20 samples, the 

cross-correlation becomes 0 as shown below. It is evident from the nature of 
trigonometric functions as orthogonal functions system. 
Zh[i]h[i-fau] (au=20 



25 =Z h [i] h [i-20] 

i=0 
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— lcos 10i + cos 2iJ + [cos 14i + cos 2i] } 

2 20 20 2 20 20 

= 0 

According to similar calculations, obviously, the cross-correlation 
1 5 becomes 0 every time when the time difference tau = ±20, ±40, ±60, ±80 as 
shown below. 

Z h [i] h p-teuj = 0 (tew= 0, ±20, ±40, ±80} 

i 

Synthesizing the foregoing results, the cross-correlation of h[i] and hp] 
20 becomes impulse form. 

40 (tau = 0) 

I h [i] h [\-tau] = 

i 0 (fau=0, ±20, ±40, ±80) 

Now, the cross-correlation of hp] and h 1 p] is determined. 
25 The value of cross-correlation becomes 10 when the time difference tau 

is 20 samples and -20 samples as shown by the following expressions. 

I hfflhp-tei/i^ 

i 
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When the time difference tau is other than ±20 samples, the cross- 
correlation becomes O from the nature of trigonometric functions as orthogonal 
20 functions system. 

The foregoing results can be synthesized as follows. 

V 9 u n : ,- , , r 10 ( tau = ±20) 
Z h [i] h [i-fau] = { — ' 

i 0 {tau= 0, ±40, ±60, ±80) 

Synthesizing all results mentioned above, the nature of the cross- 
25 correlation of h[i] and h'fj] which is h[i] right and left inversed on the time axis 
can be expressed by graph as shown in Fig. 25 and Fig. 26. 

In the foregoing graph, as continuous sweep signals hp] and h'[i] are 
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approximated by hp] and h[i] whose frequency varies gradually, the cross- 
correlation of h[i] and h[i] presents a relatively large amplitude because the 
approximation is rough; however, it can easily be understood that this value 
becomes smaller if the approximation accuracy is improved by subdividing the 
frequency variation step. 

The foregoing results show that the cross-correlation of h[i] and hp] and 
the cross-correlation of h'[ij and h*[i] are impulse, the cross-correlation of h[Q 
and h'p] becomes almost Q. 

In other words, the linear convolution of hp] and h'p] becomes an impulse, 
the linear convolution of hp] and hp] becomes substantially Q, and the linear 
convolution of h'p] and h'p] also becomes substantially 0. 

Therefore, it is an object of the present invention to provide a digital 
modem of simple communication method not requiring complicated diffusion 
symbol or cycle control taking profit of a nature that the auto-correlation and 
cross-correlation between a signal h[i] of which group delay is proportional to 
the frequency and a signal h'p] which is hp] right and left inversed on the time 
axis become impulse and substantially Q, 
[Means to Solve the Problems] 

In order to achieve such object, the present invention is composed as 
follows. 

In short, the invention of claim 1 is a digital modem, comprising : 
a modulation circuit for outputting respectively as modulation signal, 
a signal h[k] whose amplitude frequency characteristics are constant, 
and the phase thereof varies in proportion to the square of the frequency 
(namely, group delay proportional to the frequency), 
when the transmission symbol is 1, and 



a signal h[-k] which is the signal h[k] right and left inversed on the time 

axis, 

when the transmission symbol is — 1; and 

a demodulation circuit for determining the difference between, 

the square after linear convolution of this modulation signal and the 

signal h[-k] which is the sequence h[k] right and left inversed on the time axis, 

and 

the square after linear convolution of the modulation signal and the 
sequence h[k], 

for modulation / demodulation of 1 /- 1 binary signal. 

The invention of claim 2 is the digital modem of claim 1 , wherein : 

an analogue circuit is adopted for output processing of a modulation 

signal in said modulation circuit and / or convolution processing of a modulation 

signal in said demodulation circuit. 

The invention of claim 3 is the digital modem of claim 1, wherein : 
said signal h[k] is a sequence h whose discrete Fourier transform is 
DFT(h[k]) = { cos #n 2 +pn£n 2 (0^n^U2) 
cos /J (L-n) 2 - jsin j8 (l_-n) 2 (L/2<n<L) 

where, L is the length of the sequence b, the range of k is 0 < k < L, and £ is 

a constant taking a value other than 0. 

The invention of claim 4 is the digital modem of claim 1 , wherein : 
said signal h[k] is a sequence h making 

h[k] = 1 - 2mod 2 C J^!_ ) 

2L 

where, I is the length of the sequence h, the range of h is 0 ^ k < L, mod 2 (x) is 
the remainder of division of x by 2, and x ] is an integer not exceeding x. 



The invention of claim 5 is the digital modem of claiml , wherein two 
sweep signals used for modulation of the binary signal and two FIR filter 
coefficients used for demodulation are defined as 

yM = sin( a n* + 8 n + r) (0 n<N) 

y 2 [n]=y 1 [N »1 -n] (0 n<N) 

where N is length of sweep signal y^n], y 2 [n], a, £and r are arbitrary 
constants. 

The invention of claim 6 is the digital modem of claiml , wherein two 
sweep signals used for modulation of said binary signal are defined as 
equation (01) 

y t [n} = sin(an 2 + j3n + r) (0 n<N) 

y 2 [n3 = y^[N "1 -n] (0 n<N), 

two FIR filter coefficients used for demodulation thereof are defined as 
equation(02) 

hi[n] = { 1 (IfO^YiM 
"1 Ofy,[n]<0) 

h a [n] = h, [N »1 -nj (0^n<N) f 

where N is length of sweep signal y^n], y 2 [n] and FIR filter coefficient h,[n], 
h 2 [n], a , $ , r are arbitrary constants. 

[Detailed Description of Preferred Embodiments] 

Now, embodiments of the present invention will be described referring to 
drawings. 

Fig. 1 shows the signal flow of a digital modem to which the present 
invention is applied. 

The digital modem comprises a modulation circuit 1 and a demodulation 
circuit 2, for modulation / demodulation of 1 /- 1 binary signal. 
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When a binary signal of value 1 or 0 is treated, it may be simply 
considered by substituting the signal of value - 1 with the value 0. 

In the drawing, 11 of the modulation circuit 1 is a generator of sequence 
h[k] of finite length, and 12 is a generator of sequence h[-k] of finite length, 
5 which is h[k] whose time axis is inverted. 

When these circuits are realized by hardware, ROM (Read Only 
Memory) circuit may be used. 

A switch 13 is a selector for changing over h[k] and h[-k] according to 1 
or -1 input signal, and selects the output of the generator 1 1 generating h[k] 
1 0 when the input signal is 1 , and selects the output of the generator 1 2 generating 
h[-k] when the input signal is - 1 . 

At this point, by reversing the symbol, it goes without saying that h[-k] is 
generated if an input signal is 1 and h[k] is generated if an input signal is - 1 . 
Here, the generators 1 1,12 are based on digital processing and also may 
15 be analogue processing. 

In this case, for an analogue processing device, VCO and SAW oscillator 
or the like are preferably capable of generating a wide band sweep signal. 

In the first embodiment, the sequence h[k] of finite length is generated as 
follows. 



20 



exp[ji^ 1 n 2 ]^ cosily 1 rf\ + jsm[— JL n 2 ] 

La La La (O^n^L/2) 

H[n] = 



25 



H*[L-n) =cos[£^l(L-n) 2 ]-pin[^l_L L-n) 2 ] 

La La 

(L/2<n<L) 

(D 

H[n] = { cos * n * + & n * n ' (O^n^L/2) 
30 cos 0 (L- n) 2 - jsin 0 (L— n) 2 (L/2<n<L) (1 )' 
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(2) 



g[k] = reai(lDFT (H[n])) = real [. 



1 



I H[n]exp[/ 



; 2n 



kn]] (0^k<L) 



(3) 



L 



L 



g[k + 1 [1 - Jj [0^k<J^[1 - JL 
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h[k] = { 



(4) 



g[k-l [1+1]] [1[1 - J_]£k<L] 



2 a 2 a 



IDFT : inverse discrete Fourier transform 

L : penod * : complex conjugate f = - 1 

A complex sequence H[n] of length L shown in the equation (1) is inverse 
Fourier transformed as shown by the equation (3) before taking the real part to 
generate a sequence g[k], and further, g[k] is shifted circularly to the left as 
shown in the equation (4) to obtain h[k]. 

The parameter a in the equation (1), (3), (4) is a parameter to define the 
nature of h[k]. 

The value of the parameter a is enough larger than 1 , and in general, it 
may be 2. 

Examples of h[k], h[-k] generated with the length L=512 f parameter a =2 
by the equations (1), (3), (4) are shown respectively in Fig. 2 and Fig. 3. 

Here, the sequence h[k] generated by the equations (1), (3), (4) is known 
as TSP {Time Stretched Pulse) and the generation method thereof is publicly 
known. 

Besides, 21 of the demodulation circuit 2 is FIR filter having as filter 
coefficient the sequence h[-k] which is h[k] whose time axis is inverted, and 22 
is FIR filter having as filter coefficient the sequence h[k]. 

23 and 24 are square multipliers. 
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Input modulations signals are filtered and output respectively by the FIR 
filters 21, 22, squared by the square multipliers 23, 24, and the difference of 
results is determined to obtain the demodulation output. 

Here, in case where the input signal is the sequence h[k], the output from 
the FIR 21 , in short, linear convolution of h[k] and h[-k] becomes an impulse as 
shown in Fig. 4. The linear convolution of two signals is defined by the equation 
(7). 

y[n] = h[n] * x[n] 

= L Ih[k]x[n-k] (7) 

k=0 

On the other hand, the output from the FIR 22, in short, linear 
convolution of h[k] and h[k] becomes a signal of narrow amplitude as shown in 
Fig. 5. 

In case where the input signal is the sequence h[-k], the output from the 
FIR 21 becomes the signal shown in Fig. 5, and the output from the FIR 22 
becomes the signal shown in Fig. 4. 

Consequently, a pulse train having a polarity corresponding to the value 
of binary signal input of the modulation circuit 1 to output from the demodulation 
circuit 2 is obtained. 

Here, FIR filters 21, 22 are based on digital processing and also may be 
analogue processing. 

In this case, an analogue processing device has a characteristic that the 
higher the frequency is, the larger the delay time is. 

For example, CCD and SAW variance type delay lime or the like are 
preferably used. 
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In case of the first embodiment, the linear convolution of h[k] and h[-k] 
becomes an impufse, allowing to obtain the impulse response itself of the 
transfer system as demodulated output. 

Therefore, the transfer system characteristics can be compensated using 
5 the demodulated output during signal reception. 

In the second embodiment, a sequence h[k] generated by the expression 

[5] is used. The sequence h[k] generated here corresponds to the h[k] 

generated by the equations (1), (2), (4) , clipped and two-valued, and has the 

nature of cross-correlation and linear convolution similar to the h[k] generated 

10 by the equations (1), (2), (4) . 

h[k] = 1 - 2mod 2 [ ] {5) 

2L 

mod 2 (x) : x modulo 2 
: floor 
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x-mod 2 (x) =n (6) 



2 

n : integer 



in the equation (5), mod 2 represents the remainder operation taking 2 as divisor. 
20 x is an integer not exceeding x, namely operator for rounding the 

real x to an integer. 

The length of the sequence to be generated is specified by L 
The generated signal becomes a binary value having the value 1 or -1 . 
Examples of h[k] generated by the equation (5), taking 40 as the value of 
25 L and h[-k] which is h[k] whose time axis is inverted are shown respectively in 

Fig. 6 and Fig. 7. 
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Here, in case where the input signal of the demodulation circuit 2 is the 
sequence h[k] } the output from the FIR 21, in short, linear convolution of h[k] 
and h[-k] becomes a signal in impulse form as shown in Fig. 8. 

On the other hand, the output from the FIR 22, in short, linear 
convolution of h[k] and h[kj becomes a signal of small amplitude as shown in 
Fig. 9. 

In case where the input signal is the sequence h[-k], the output from the 
FIR 21 becomes the signal shown in Fig. 9, and the output from the FIR 22 
becomes the signal shown in Fig. 8. 

In case of the second embodiment, h[k] and h[-k] being binary signal, it 
can be realized by an extremely simple circuit, in case of hardware embodiment. 

To be more specific, examples of modulation signal and demodulated 
output signal will be shown below. 

In case where the input signal of the modulation signal 1 is a 8 bit {1,-1, 
1,1,-1,1,-1,-1} shown in Fig. 10, the output from the modulation circuit 1 will 
be as shown in Fig. 1 1 . 

When the modulation signal shown in Fig. 1 1 is input to the demodulation 
circuit 2, the demodulated output becomes a signal in pulse form shown in Fig, 
12. 

The pulse polarity of demodulation output corresponds to the value of the 
input signal of Fig. 10. 

A bit train of value 1 or -1 can be restored easily from the signal of Fig. 
12, by employing a demodulation circuit 2 using DLL (Delay Locked Loop). 

In the third embodiment, two sweep signals used for modulation are 
generated on a basis of the definition of equation(01). 
y,[n] = sin(an 2 + £n + r) (0^n<N) 
YzM = ft [N - 1 -n] (0^n<N) 
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Examples of the two sweep signals y^n], y 2 [n] generated by designating 
values of a, j3, r as 0.001, 0, 0 in equation(01) are respectively shown in 
Fig.28, 29. 

Further, the result of linear convolution processing of y,[n], y 2 [n] are 
shown in Fig.30. The result of linear convolution processing of y^n], yjn] are 
shown in Fig,31. 

From the result of linear convolution of Fig.30 and Fig.31, if y^n] and, 
y 2 [n] are used for modulation and demodulation, the circuit using FIR filter taken 
VM 5 y 2 [n] as coefficients is capable of demodulating a signal. 

The results of modulation / demodulation are shown in Fig.36 and Fig.37. 

Fig,36 is a modulation waveform using y,[n] , y 2 [n] and Fig.37 is a 
demodulation waveform thereof. 

In the fourth embodiment, two sweep signals y^n] , y 2 [n] used for 
modulation are generated on a basis of the definition of equation(01). 
y^n] = sin( a n 2 + jS n + r) (0^n<N) 
y 2 [n]=y, [N -1 -n] (0^n<N) 

Two FIR filter coefficients h^n] , h 2 [n] used for demodulation are defined 
as equation(02). 

hiM = { 1 Ofo^M 
-1 (ifVi[n]<0) 
h 2 [n] = h 1 [N -1 -n] <0^n<N) , 

Examples of two sweep signals y t [n] , y 2 [n] generated by designated 
values of a, £ , r as 0.001,0,0 in equation(01) are shown in Fig.28, Fig.29. 

Further, FIR filter coefficients h t [n] , h 2 [n] generated on a basis of the 
definition of equation(02) are shown in Fig.32, Fig.33. 

Furthermore, the results of linear convolution processing of yjn] and 
hjn] are shown in Fig.34. The results of linear convolution processing of y^n] 
and hi[n] are shown in Fig.35. 
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From the results of linear convolution processing of Fig.34 and Fig.35, if y n [nl 
and y 2 [n] are used for modulation, the circuit using FIR filter taken y t [n] and 
y 2 [n] as coefficients is capable of demodulating a signal. The results of 
modulation / demodulation are shown in Ffg.36 and Fig.38. Fig.36 shows a 
modulation waveform using yjn] and y 2 [n]. Fig.38 shows waveform of a signal 
demodulated using FIR filter taken hjn] and h 2 [n] as coefficients. 

As mentioned above, the digital modem of the present invention 
comprises a modulation circuit for outputting as modulation signal either the 
signal h{k} whose amplitude frequency characteristics are constant, and the 
phase thereof varies in proportion to the square of the frequency (namely, group 
delay proportional to the frequency) or the signal h[-k] which is the signal h[k] 
right and left inversed on the time axis according to a binary signal of 1 /- 1 and 
a demodulation circuit for determining the difference between the square after 
linear convolution of this modulation signal and the signal h[-k] which is the 
sequence h[k] right and left inversed on the time axis, and the square after 
linear convolution of the modulation signal and the sequence h[k], outputting as 
a demodulation signal an impulse signal train having the puise polarity 
corresponding to the binary signal of the transmission signal. 

Consequently, the present invention allows to realize a digital processing 
using the sequence h[k] whose characteristics are strictly defined and, at the 
same time, to simplify the apparatus and thereby improve the reliability and 
economy, because complicated generation circuits of diffusion symbol required 
conventionally for the transmitter and the receiver become both unnecessary. 

In addition, the demodulation circuit 2 can be simply composed only of 
digital filters such as DSP having the sequence h[kj as filter coefficient, 
complicated synchronous control for receiving a signal being unnecessary. 
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Further, it is resistant to intermittent pulse noise, as the modulation / 
demodulation is performed by a kind of pulse signal expansion/compression on 
the time axis. 

Such resistance to the noise is described below. 

As shown in Fig. 27, suppose that input signal of the demodulation circuit 
y(ij, outputs from two FIR digital filters respectively are p1[i], p2[i], outputs from 
two square circuits respectively p^fij, p 2 2 p], and output from the demodulation 
circuits wfj]. 

Here, as shown by the following equation, the input y[i] is supposed to be 
the sum of the modulation signal hp] and additional random noise np]. 
y[i] = h(i] + n{i] 

P1PL Pi 2 [i] are determined as follows. Here, * is an operator representing 
the linear convolution, h'[ij is sequence h[i] whose time axis is inverted, and 6 
[i] a digital function. 

As described already, the linear convolution of hp] and h'p] becomes 
substantially a delta function. 

PiP] = yp]*h'[i] = (hp] + n[i])* h'[i] 

= h[i]*h'P] + np]*h'fi] 
= d p] + np] * h'p] 
Pi 2 [i] = (<5[i] + n[i]*h'p]) 2 

= S 2 p] + 2 «5 [i] • n[Q * h'[i] + (np] * h'ffl) 2 
Next, P 2[i], p22[i] are determined as follows. Here, m[i] is the linear 
convolution of hp] and hp] and, as described already, is a low power signal 
whose amplitude is near 0. 

Psffl = y[i] * hp] = (hp] + np] ) * hp] 
= h[i]*h[i] + n[ij*h[i] 
= mp] + np] * hp] 
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p 2 2 M = (m[i] + n[i]*h[iJ) 2 

= m 2 [i] + 2m[i] • np] * hp] + (n[0 * h[i]) 2 
Next, the output w[i] of the demodulation circuit is determined as follows. 
w[ij = p 1 2 [i] -p 2 2 [ij 

= 6 2 [i] + 2 6 p] • n p] * h'[i] + (n[i] * h'[i]) 2 ~m 2 [i] - 2m[i] ■ n[i] * h[i] - 
(n[i]*h[i]) 2 

= <5 [i] + 2 <5 [i] ■ n[i] * h'[ij -m 2 [0 - 2m[i] • nfj] * h[i] + (n[i] * h'p]) 2 - 
(n[i]*hpj) 2 

= <5 p] (1+ 2 ■ np] * h'pj) -m 2 p] - 2mp] ■ np] * hp] + (np] * h'p]) 2 - 

If the random noise np] is stationary in an interval shorter than tine length 
of hp] and h'p], the squares of the linear convolution of hp] and h'p] become 
equal each other. 

(n[i]*hp]) 2 -(n[i]*h'P]) 2 
Consequently, w[i] becomes as follows, 
wpj * d [ii (1 + 2 • nffl * h'p]) -m 2 [i] - 2mp] • np] * h[i] 
The value of the second term and the third term of the light side of the 
foregoing expression, being the product of narrow amplitude signals each other, 
can be taken as substantially 0. 
|m[i]|< 1 
[n[i] * h'PJl < 1 
Hi]*h[i]|<1 
Consequently, wp] becomes further as follows. 
w[i]* «Jp](H2 • n[0 * h'[i]) 

Here, k, is a smalt constant whose absolute value is equal or inferior to 1. 

IM<1 
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Therefore, it is understood that the modulation / demodulation method of 
this method can perform the demodulation processing without being 
substantially influenced by noise when random noise is applied to the transfer 
path, and only the pulse amplitude of the demodulation output varies. 

Besides, in case where the input of the demodulation circuit is the sum of 
the modulation signal h'[i], namely hp] whose time axis is Inverted and the 
random noise n[i], w[i] becomes as follows. 

In this case also, it is understood that only the negative pulse amplitude 

varies. 

w[i] - 6 [i] (1 - 2 • n[i] * h[i]) 

= -<5[i]{1+k 2 ) 
N<1 

Further, Fig. 13 to Fig. 20 shows examples of demodulation output 
waveform, in case where random noise is actually applied to the modulation 
wave in the transfer path. 

BRIEF DESCRIPTION OF THE SEVERAL VIEWS OF THE DRAWING 

Fig. 1 shows the signal flow of a digital modem to which the present 
invention is applied; 

Fig. 2 shows h[k] generated by equations (1), (3), (4) supposing L=512 
and a =2; 

Fig. 3 shows a sequence h[-k] which is h[k] shown in Fig. 2 whose time 
axis is inverted; 

Fig. 4 are linear convolution of h[k] shown in Fig. 2 and h[-k] shown in Fig. 

3; 

Fig. 5 are linear convolution of h[k] and n[k] itself shown in Fig. 2, or 
linear convolution of h[-k] and h[-k] itself shown in Fig. 3; 

Fig. 6 is h[k] generated by the equation (5) supposing L=40; 
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Fig. 7 shows a sequence h[-k] which is h[k] shown in Fig. 6 whose time 
axis is inverted; 

Fig. 8 are h[k] shown in Fig. 6 and linear convolution of h[-k] shown in Fig. 

7; 

5 Fig. 9 are h[k] and linear convolution of n[k] itself shown in Fig. 6, or h[-k] 

and linear convolution of h[-k] itself shown in Fig, 7; 

Fig. 10 is an example of input signal to the modulation circuit, 8bit serial 
signal of {1,-1, 1,1, -1,1, -1,-1}; 

Fig. 11 shows a modulation output signal when the signal of Fig. 10 is 
10 applied as input to the modulation circuit of Fig. 1; however, here, the 
generators 11, 12 are supposed to output the binary signals shown respectively 
in Fig. 6, Fig. 7; 

Fig. 12 shows the demodulation output signal when the modulation 
signal of Fig. 11 is applied as input to the demodulation circuit of Fig. 1; 
15 however, here, the coefficients of FIR filters 21, 22 are supposed to be binary 
values h[-k], h[k] shown respectively in Fig. 7, Fig. 6; 

Fig. 13 is an example of modulation waveform of the first embodiment, 
wherein the modulation is performed by the circuit of Fig. 1 , by taking as input a 
9bit serial signal #{-1,1,-1,-1,1,1,1,-1,1}. 
20 The sequence h in Fig. 1 uses the one generated by expressions (1), (3), 

(4). 

Here, parameters in the expressions are set to L=64 and ct =1 ; 

Fig. 14 shows a demodulation output when the signal of Fig. 13 is input 
to the demodulation circuit; it is understood that a pulse train corresponding to 
>5 the bit train {-1,1,-1,-1,1,1,1,-1,1} input to the modulation circuit of Fig. 1 is 
output; 

Fig. 15 shows a waveform on evenly distributed noise; 
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Fig. 16 shows a signal waveform wherein noise of Fig. 15 is added to the 
modulation waveform of Fig. 13; 

Fig. 17 shows the demodulation output of the signal of Fig. 16 by the 
circuit of Fig. 1 , demonstrating that the modulation / demodulation method of the 
present invention is resistant to the noise, because the demodulation 
processing is performed normally even when random noise is added to the 
signal; 

Fig. 18 shows a signal wherein sinusoidal waves of cycle 5.9 sample and 
cycle 13.2 sample are added; 

Fig. 19 shows a signal waveform wherein the signal of Fig. 18 is added 
to the modulation waveform of Fig. 13; 

Fig. 20 shows the demodulation output of the signal of Fig. 19 by the 
circuit of Fig. 1 , demonstrating that the modulation / demodulation method of the 
present invention is resistant to the parasite, because the demodulation 
processing is performed normally even when cyclic parasite signal is applied; 

Fig. 21 shows a signal waveform of the sweep signal h[i]; 

Fig, 22 shows a signal waveform of the sweep signal h'p]; 

Fig. 23 shows a signal waveform approximating the sweep signal h[i] by 
a cos signal whose frequency varies gradually; 

Fig. 24 shows a signal waveform approximating the sweep signal h'[i] by 
a cosine signal whose frequency varies gradually; 

Fig. 25 

is a graph showing that the nature of the cross-relation of hp] and h[t] and 
h'p] and h'[i] takes an impulse form; 
Fig. 26 

is a graph showing that the nature of the cross-relation of hp] and h'{i] 
becomes substantially 0; and 
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Fig. 27 is signal processing diagram illustrating that the demodulation 
circuit of the present invention is resistant to noise. 

Fig. 28 is an example of sweep signal y,[n] generated by designating 
respectively values of a , £ and r as 0.001 , 0 and 0 in equation(01 ). 

Fig.29 is an example of sweep signal y 2 [n] generated by designating 
respectively values of a , £ and r as 0.001 , 0 and 0 in equation(01). 

Fig.30 shows linear convolution of y,[n] shown in Fig.01 and y 2 [n] shown 
in Fig.02. 

Fig.31 shows linear convolution of y^n] shown in Fig.01 and y^n] itself. 
Fig.32 shows FIR filter coefficient h,[n] generated on a basis of the 
definition of equation (02). 

Fig.33 shows FIR filter coefficient h 2 [n] generated on a basis of the 
definition of equation{02). 

Fig.34 shows linear convolution of y,[n] shown in Fig.01 and h 2 [n] shown 
in Fig.06. 

Fig.35 shows linear convolution of y,[n] shown in Fig.01 and h,[n] shown 
in Fig.05. 

Fig.36 shows modulation waveform using sweep signals y,[n] and y 2 [n] 
generated on a basis of the definition of equation(01). 

Fig.37 shows waveform of a signal demodulated using FIR filter waveform 
y,[n] and y 2 [n] generated on a basis of the definition of equation(01). 

Fig.38 shows waveform of a signal demodulated using FIR filter 
coefficients h^n] and h 2 [n] generated on a basis of the definition of 
equation(02). 
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